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Summary. — For stationary vacuum spacetimes the Bianchi identities of the second 
kind equate the Simon tensor to the Simon-Mars tensor, the latter having a clear 
geometrical interpretation. The equivalence of these two tensors is broken in the 
nonvacuum case by additional source energy-momentum terms, but absorbing these 
source terms into a redefinition of the Simon tensor restores the equality. Explicit 
examples are discussed for electrovacuum and rigidly rotating matter fields. 



1. Introduction 

Using the 1 + 3 splitting of spacetime associated with an observer congruence u (time- 
like, unit tangent vector to the congruence), the Bianchi identities for the curvature tensor 
can be cast into a Maxwell-like form for the electric and magnetic parts of the Weyl ten- 
sor with sources arising from the Ricci tensor terms, in turn re-expressed in terms of the 
energy-momentum tensor of the spacetime through the Einstein field equations P3 El E| • 

For the stationary case with u directed along the timelike Killing vector field, these 
equations are greatly simplified and suggested the introduction of the (complex) Simon 
0] and Simon-Mars tensors |SJ [7| EJ E3 \U\ ■ The Simon tensor was introduced 
for vacuum stationary spacetimes as a complex generalization of the Cotton- York tensor 
1 121 1 1 3 j . to which it reduces in the static case. The Simon-Mars tensor, instead, is 
a natural algebraic combination of the self-dual fields associated with the Papapetrou 
field of the timelike Killing vector field and the Weyl tensor. In the vacuum case the 
Bianchi identities state the equivalence between these two tensors, and the vanishing of 
the Simon-Mars tensor has been understood as corresponding to the alignment of the 
principal null directions of the Weyl tensor of the spacetime with those of the Papapetrou 
field associated with the timelike Killing congruence. This is exactly what happens in 
the case of the Kerr solution, uniquely characterized by the vanishing of the Simon (and 
Simon-Mars) tensor. 
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In the nonvacuum case the equivalence between the Simon and the Simon-Mars tensor 
is broken because of the presence of Ricci tensor terms in the Bianchi identities which 
may be re-expressed in terms of the source energy-momentum tensor by the Einstein 
equations. To characterize exact solutions in a way similar way to the vacuum case, one 
approach would be to study these distinct tensors to explore their meaning. Alterna- 
tively, one can absorb the source terms into a redefinition of the Simon tensor so that 
the resulting generalized Simon tensor equals the Simon-Mars tensor, which continues 
to have a clear geometrical meaning. In this more geometrically motivated case, one 
can classify exact solutions studying only the vanishing of a single tensor. Here we ap- 
ply these considerations to electrovac stationary spacetimes which contain a source-free 
electromagnetic field ^T] and stationary spacetimes containing a rigidly rotating perfect 

fluidEH. 

2. The 1 + 3 Maxwell-like equations for the Weyl tensor with sources 

The Bianchi identities of the second kind equate the divergence of the Weyl curvature 
tensor C a p 7 s to half the Cotton tensor which can be thought of as a current in 
analogy with electromagnetism ^5] (where F a ^.^ = 4irJ a ) 

(1) c aS p~ n5 = -v [/3 (Wir 7] - i (4) i?^ 7] ) = J a 01 , 

following the conventions of Misner, Thorne and Wheeler including the signature +2. 
It is convenient to single out the Ricci tensor (RT) and the Ricci scalar (RS) contributions 
to the Cotton current J a p-y of Eq. QJ: 

(2) J a 'pj = J(nT) a pj + J(RS) a i3-r, 
where 

( 3 ) ^(RT) Q /3 7 = - V [/3( (4) -R Q 7 ]), ^(RS) Q /3 7 = g V [/3 (4) ^ tt 7]- 

The Einstein equations (in Ricci form) 

(4) R a p = n(T a p - ±T\5 a p) 

can then be used to replace the Ricci tensor and scalar curvature terms in this current 
by the energy-momentum tensor: 

(5) J a Pl = - K V [/3 (T Q 7] - \t6%$. 

Splitting these equations with respect to a generic time-like congruence u leads to their 1 + 
3 Maxwell-like form 3] given explicitly in our notation in ^3 . The Weyl tensor splits 
into two symmetric tracefree spatial fields, its electric and magnetic parts respectively 



(6) 



E a p = C a 7 p S u^u 5 , H a p = -*C%p S W<u s = ^ n V CPsfip u p 
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If a is the acceleration vector and uj the vorticity vector of u, introducing the complex 
spatial fields from the gravitoelectromagnetic connection and Weyl curvature fields 

(7) z = -a-iw = g-iH/2 , Z = E - iH , 

where g and H are the gravitoelectric and gravitomagnetic vector helds, leads to a 
more efficient representation of those equations. Here we specialize to the case of a 
stationary spacetime and align u = A/ _1 £ with the associated Killing vector £, so that 
the acceleration can be expressed as a = V InM, where M = \ — S, a S, a \ 1 ^ 2 - 

Splitting the Weyl divergence equations and expressing them in terms of an adapted 
frame {ey = u, e a }, a = 1, 2, 3, where the spatial frame {e Q } is orthogonal to u, one finds 
in an index-free notation the identities |10| 

(8) Af 3 div[Af- 3 Z] = 3z ■ Z + p {G) , NT 1 Scurl[MZ] = z x Z + i.J {G) , 
where the complex current fields are defined by 

(9) p{ G) = J T aT + U* T aT , = -[J(ab)T + U\ab)T] ■ 

The div and Scurl spatial operators on symmetric spatial 2-tensors S 

(10) [div(u)S*] Q = V(u),gS Q/3 , [Scurl(w)5] Q/3 = r,(u)*V (u) 7 5 /3) s 

are defined in terms of the spatial covariant derivative V(ii), the spatial volume 3-form 
Vi^afl-y — vfrfsapj (associated with the spatial duality operation *< u >, used for example 
to define the vorticity vector ui(u) a = ^ry(u) a ' 37 V i au 7 from the corresponding 2-form) 
and the spatial cross and dot products associated with u, all given explicitly in jl()| by 
spatially projecting spacetime derivatives and quantities orthogonally to u. For example, 
the cross-product of a spatial vector and spatial symmetric tensor is 

(11) [X x u A] a(} =i 1 (u)^ a X 1 A^ s ■ 

For simplicity we drop the reference to u on these spatially projected operators, resorting 
to the abbreviation V = V(u) to distinguish the spatial covariant derivative from the 
spacetime covariant derivative V, an unnecessary distinction when acting on stationary 
scalars. In an adapted frame spatial quantities may be expressed using only Latin indexed 
components. 

In the Scurl equation in the projected Bianchi identity JSJl we have three terms: 

1 . The (rescaled) Scurl of Z which defines the complex, symmetric and tracefree Simon 
2-tensor ^ of the vacuum spacetime 

(12) Simon(vacuum) = AT 1 Scvx\[MZ]. 

2. The spatial cross product of z and Z which defines the complex, symmetric and 
tracefree Simon-Mars 2-tensor [3 El HI HO] 



(13) 



SimonMars = z x Z 
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3. The Cotton gravitational current term iJ^ G \ which is also complex, symmetric and 
tracefree. Its expression in terms of Ricci tensor components is 

= ^(aM 2 ^R b)T \{ U x ( 4 >Ricci) Qb , 

(14) ~(Scurl( 4 >Ricci) Qb + |[ W(a ( 4 )i? b)T ]( TF ) . 

Here ( 4 ) Ricci denotes the index- free form of ^'R a b to distinguish it from the spacetime 
curvature scalar ( 4 ' R. The only nonvanishing field arising from splitting the Ricci scalar 
contribution to the Cotton current is the spatial tensor 

(15) J {KS) a bc = l -d [b ^R5 a c] , 
so that 

J(RS) a fcT = , *J (RS) V - —rffdJVR . 

Thus * J(Rs)(ab)T = and hence it does not contribute at all to J^ G \ 

Let us define the nonvacuum Simon 2-tensor for any stationary spacetime (with any 
energy-momentum source) as 

(16) Simon(nonvacuum) = Simon(vacuum) — i J^ G ' , 

so that the Simon (nonvacuum), hereafter denoted only by Simon, and the Simon-Mars 
tensors also remain equal in the nonvacuum case. Of course once the Ricci tensor in the 
source term is replaced using the Einstein equations (J3J , this equality between Simon 
and Simon-Mars is only true if they are each evaluated on a solution of those equations, 
in contrast to the vacuum case where they are identically equal for any metric. With 
this understanding there is no longer any need to use two different names, Simon and 
Simon-Mars, to denote the same quantity. 

It is worth noting further that this identification has a true geometrical content. In 
fact, the vanishing of the Simon-Mars tensor (z x Z — 0) is equivalent to the alignment 
of z and Z, in the sense that Z must be proportional to the tracefree tensor product of 
z with itself [Ej 

(17) Zoc [z(g>z]( TF) . 

Of course, other generalizations of Simon(vacuum) can also be studied. 

Finally, we remark that when working with explicit solutions of the Einstein field 
equations, 

J(G) 

can have useful special forms. For example, when = — iVlna x Z, 
then the Scurl identity for any scalar a and any symmetric spatial 2-tensor A 

(18) aScml[a- 1 A] = -Vina x A + Scurl A 
enables this current to be incorporated into the existing Scurl term 



(19) 



Simon = (cr/M) Scurl[o- _1 MZ] 
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This again would re-establish the vacuum property that the Simon tensor is proportional 
to a Scurl. This is the case with the family of Kerr-Newman-Taub-NUT electrovac 
spacetimes (including the more physically interesting Kerr-Newman spacetime), where 
the extra factor of a converts this to the same Scurl expression when written in terms 
of the Ernst potential as in the vacuum formula (modulo a scale factor). For rigidly 
rotating perfect fluid solutions, as in the case of the Kramer solution |18|. is exactly 
zero and the vacuum case formulas again hold. 

3. Electrovac spacetimes: Maxwell equations and the Papapetrou fields 

The stationary Maxwell-like equations © for the Weyl tensor expressed in complex 
self-dual form are very similar to the stationary Maxwell equations expressed in complex 
form ^D] 

(20) M~ 2 div(M 2 x) = -2z-x + 4tt P , NT 1 curl(Mx) = 4niJ 

satisfied by the spacetimc's electromagnetic field (using an arrow to distinguish the elec- 
tric field from the electric part of the Weyl tensor in index-free notation) 

(21) F = uAE + *™B , x = E-iB. 

They only differ in form by the missing z x x term in the curl equation here, where *< u iB 
is the spatial dual of the magnetic 1-form field B. However, when z and Z are aligned 
(i.e., zxZ = 0), the correspondence is complete. The sources p and J for the spacetime 
electromagnetic field are zero for the electrovac spacetimes considered here. The curl 
equation for the electromagnetic field implies that Mx is the gradient of a potential 

(22) Mx = V$ , or x = M~ X V$ . 

Papapetrou ^H] observed that a Killing vector field £ can be interpreted as a Lorentz 
gauge vector potential for an electromagnetic field later called the Papapetrou field [HI E] 

(23) F af3 = [d£] a/3 = - Vp£ a = 2V Q £ /3 

satisfying Maxwell's equations with an effective current source proportional to the Ricci 
tensor 

(24) £ (a;/3) = - £% = - Ca;P' P - -&R? a ~> £»/J ! " = ^RPc , 

where this divergence condition is merely the trace of the Ricci identity for a Killing 
vector field. Rescaling the Papapetrou field F by the factor (2M) _1 ^U] and introducing 
its complex self-dual representation yields the self-dual 2-form T associated with z 

(25) T = (2M)-\F + i*F) = u A z + i*M Z , 

which when multiplied by M satisfies the same Maxwell equations but with nonzero 
sources generated by the electromagnetic field. These equations 



(26) 



M~ 3 div(M 3 z) = -2z ■ z - x ■ x , M~ 2 curl(M 2 z) = x x x 
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are obtained by splitting the vector Eqs. (|24|l and using the Einstein equations (J2J to 
replace the Ricci tensor appearing there in terms of the appropriate projections of the 
electromagnetic energy-momentum tensor of the spacetime 

(27) kT T j- = —x ■ x , inT T a — [x x x\ a , 

which then leads to the source terms appearing in these Maxwell equations. 

Inserting the potential representation (IL'l'l) of Mx into the rescaled Papapetrou field 
curl equation of Eqs. l|2"6*)l then yields 

(28) M~ 2 curl(M 2 z) = M" 2 V$ x V$ = -A/" 2 curl($V$) , 
namely 

(29) curl(Af 2 z + l>V$) = . 

This latter quantity in parentheses must therefore be a gradient, defining the Ernst po- 
tential S [T7j (see their Eq. (18.34) with F = -M 2 , reducing to the integrated Eq. (18.29) 
in the vacuum case $ = 0) modulo a numerical factor 

(30) VS = -2{M 2 z + i>V$) . 

The remaining divergence equation of Eqs. I|26l) is the Ernst equation, which can be 
expressed in the form (equivalent to Eq. (18.38) of (TJ|, see Appendix A) 

(31) Mdiv(MV£) = -2Af 2 z ■ V£ . 

This result follows from using the definition of £ to re-express x in terms of the gradient 
of $ in the two divergence equations 1)20(1 and (|2fi(l 

(32) idivV£ + #divV$ = A/ 2 (a + 2z)-z, div V$ = -(a + 2z) ■ V$ 

and using the second of these equations to eliminate div V$ in the first, yielding 

(33) divV£ + a- V£ = -2z-V£ , 

which can then be rewritten as Eq. I|31|l . 

The Simon-Mars tensor can be evaluated in a form which shows its relation to the 
Ernst potential expression for the Simon tensor, assuming that the Einstein field equa- 
tions are satisfied so that the two tensors are in fact equal. Following the notation of 
|l()j . the purely spatial components and the mixed time-space components of the Ricci 
tensor field equations (see J7|, Eqs. (18.15) and (18.16)) are respectively 

[( 4 )Ricci]iI F )=-2Af- 2 [V (a $V &) $]( TF ), 

(34) (4) i? T & = iAT 2 [curl(A/ 2 z)] b = -iAf- 2 [curl(§V$)] fe , 
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where z and $ are related to the Ernst potential by l|30|) . A straightforward computation 
gives the electric and magnetic part of the Weyl tensor in terms of the gravitoelectro- 
magnetic fields a and ui and the Ricci tensor of the quotient 3-manifold 

E = --[4w<g> 10- Va - a® a- (3) Ricci] (TF) , 

(35) H = - SYM[Vw + 2a <g> w] (TF) . 

Then using the formula 

(36) [ (3) Ricci] < TF > = [<*> Ricci + Vfl + fl ® fl + 2 w ® w ] < TF > , 
one finds 

(37) Z = E — iH =[- SYM(Vz) + z <g> z + i (4) Ricci] (TF) . 
The Simon-Mars tensor is then 

(38) SimonMars = z x (- SYM Vz + i (4) Ricci) 

since z x (z <g) z) is identically zero. By next replacing the first term here using the Scurl 
identity valid for an arbitrary spatial vector field X 

(39) Scurl(Y ®I) = -Ix SYM(VY) + ^[SYM(Y <g> curlY)] (TF) , 

and using the Einstein Eqs. I|34|) to replace the spatial components of the spacetime Ricci 
tensor, the Simon-Mars tensor becomes 

SimonMars = Scurl(z ® z) - -[SYM (2 ® curl z)] (TF) 

(40) -M~ 2 z x SYM(V$ (g) V#) . 

Finally by substituting z = M~ 2 (M 2 z) into the curl term here expanding the result with 
the appropriate product rule and incorporating the resulting additional acceleration term 
into the Scurl term, and then using Eq. (|3Tj|) 

(41) curl(M 2 z) = - curl($V$) 
to replace the other term in this expansion, one finds 

SimonMars = M~ 3 Scurl(M 3 z (g> z) - -[SYM(z <g> (x x x)] (TF) 

(42) -z x SYM(x ® x) . 

In the vacuum case where x — 0, this reduces to 

(43) SimonMars = A/~ 3 Scurl(M 3 z ® z) , 
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which can be rewritten in terms of the Ernst potential by making the substitution z = 

-\M- 2 V£ 

(44) SimonMars = ^Ar 3 Scurl(M -1 V£ ® V£) . 

As shown in , this is proportional to the form of the Simon tensor given in |17j and 
to the original form given by Simon. 

4. Rigidly rotating stationary fluids 

Consider a perfect fluid at rest in the reference frame adapted to the timelike Killing 
direction u. The energy- momentum tensore can be written as 

(45) = (p + p)u a uP +pg a P 
or in components 

(46) Ttt = p , Tm = , T ab = phab , 

where the trace is T = T a a = ip — p. The conservation equation V ' pT a P = gives 

(47) (p + p)VlnM + Vp = , 
or assuming p = p(f) with / = M 2 

< 48 » ^+§H- 

In this case one can immediately verify that = 0. Indeed the result follows from 
Eq. (fHEjl using the fact that Ricci a b is a pure trace term (killed by exterior product and 
Scurl operations) and R-j a vanishes identically. 

5. Explicit examples 

The Kerr- Newman- Taub- NUT spacetime 

In Boyer-Lindquist coordinates (x a ) — (t,r,9,<f>), with the abbreviations (c, s) — 
(cos 9, sin 9), the (exterior) Kerr-Newman-Taub-NUT (KNTN) spacetime line element is 

m 

ds 2 = --(A - a 2 s 2 )dt 2 + -[A X - a(E + a X )s 2 ]dtd(f> 

(49) +![(£ + a X fs 2 - x 2 A]d0 2 + ^-dr 2 + Ed9 2 , 

and the corresponding electromagnetic Faraday tensor can be expressed in terms of the 
2-form 

F = $r{[r 2 -(£+ ac) 2 ]dr A (di - xM) 

(50) +2rs(£ + ac)d9 A [(r 2 + a 2 + l 2 )d<f) - adt}} , 
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which is the exterior derivative of the vector potential 

(51) A = -^(d<- X d^) . 
Here S, A, and \ arc defined by 

(52) £ = r 2 + (I + ac) 2 , A = r 2 - 2Mr - I 2 + a 2 + Q 2 , X = as 2 - 2£c . 

Units are chosen so that G = c = 1, so the parameters (Ai, Q, a, i) all have the dimension 
of length. The source has mass Ai, electric charge Q, angular momentum J = Aia 
(i.e., gravitomagnetic dipole moment) along the 8 = direction, and gravitomagnetic 
monopole moment p = —I, where t is the NUT parameter. 

The unit timelike vector aligned with the Killing direction u = -rjdt (M~ 2 = S/(A — 
a 2 s 2 ) has nonzero acceleration and vorticity (if a / 0), summarized by the complex 
vector 

,„v {[r-i(i + ac)]{M+il)-Q 2 } 

(53) 2 = [r + i (l + ac)] 2 (A-a 2 S 2 )[r- i (l + ac)] { - Adr + ^ ' 
while the electric and magnetic parts of the Weyl tensor are summarized by 

< 54 > z =w^wBmm^^^ m - 

so that the Simon-Mars tensor (and hence the Simon tensor as well) vanishes. 
The Kramer solution 

Using a system of coordinates (t, /, f2, </>), directly related to the Ernst function £ = 
f + if2, the rigidly rotating perfect fluid solution due to Kramer (see Eq. 25 of JB|> with 
c = 1) is given by 

2 «j, 1 j,n2 , 1 ( 2 .,2 , 2 Jri 2 , P - Fq 2 



(55) ^ = - m - + 7 ^_d/- + ^d^ + ^d0 

where f3 = 0(f) = — i*b +d(e^//) and d, Fq are constants. The energy-momentum tensor 
of the fluid is 

(56) T=[p(f)+p(f)]u®u + p(f)g, u =^f dt 
and the energy density and the pressure depend only on / 

( 57 ) P=§(3/- 2 ), P=§(2"/) 
and are related by the equation of state 

(58) p + -ip = F //s, 
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where k is the Einstein equation constant. The four velocity of the fluid u is aligned with 
the timelike Killing direction and has nonzero acceleration and vorticity, so that 

(59) z = — a — iuj = — -^(df - idui). 

The electric and magnetic parts of the Weyl tensor are summarized by 



4/ 



1 ?,def — 4Fnf 

3 del _ 2F / d/ ^ d/ ~ ' (d/ ^ ^ + ^ ® d/) 



-^rj(3de f - 2F f)dQ ® dft + f^de' - 2F /)dw ® dw 
(60) = /[z®z]< TF >. 

A straightforward calculation shows that the Simon-Mars (and hence the Simon tensor 
too) is identically zero. The Cotton current is also zero as shown above. 

6. Conclusions 

The gravitoelectric description of the Simon and Simon-Mars tensors which play a 
crucial role in characterizing the Kerr spacetime has been generalized to the class of 
stationary electrovacuum spacetimes and stationary rigidly rotating perfect fluid sources. 
Again the alignment of the principal null directions of the Weyl and the Papapetrou is 
equivalent to the vanishing of the Simon-Mars tensor, which in turn forces the generalized 
Simon tensor to vanish as well. Examples concerning the electrovac Kerr-Newman-Taub- 
NUT spacetime and the perfect fluid Kramer spacetime have been presented. 
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